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Abstract 

In a hard- wall model of holographic QCD, we find that nonlinear boundary dynamics are required 
in order to maintain the correct pattern of explicit and spontaneous chiral symmetry breaking 
beyond leading order in the pion fields. With the help of a field redefinition, we relate the requisite 
nonlinear boundary conditions to a standard Sturm-Liouville system. Observables insensitive to 
the chiral limit receive only small corrections in the improved description, and classical calculations 
in the hard-wall model remain surprisingly accurate. 
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I. INTRODUCTION 



Holographic QCD is an extra-dimensional approach to modeling hadronic physics at 
low energies l|-l6|. Hadronic resonances are interpreted as Kaluza-Klein modes of five- 
dimensional (5D) fields with quantum numbers of the corresponding QCD states. Motivated 
by the AdS/CFT correspondence THSj], in hard- wall models 2tW\ the background is chosen 
to be a slice of Anti-de Sitter space (AdSs), with metric 



where rj^^i, is the Minkowski metric with mostly negative signature in 3-1-1 dimensions, R is 
the AdS radius, and e provides a short-distance cutoff in the model. The size of the extra 
dimension depends on Zm-, which sets the Kaluza-Klein scale and is related holographically 
to the confining scale of QCD. Alternative spacetime backgrounds have been motivated by 
D-brane configurations in string theory which give rise to QCD-like theories with chiral 
symmetry breaking and confinement, as in the Sakai-Sugimoto model based on the D4-D8 
system 10|. Non-normalizable {i.e. infinite-action) backgrounds of fields act as sources of 
corresponding operators in QCD, and normalizable {i.e. finite-action) backgrounds deter- 
mine the expectation values of those operators 
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12l |. The established AdS/CFT dictionary 



between physics in 3-1-1 and 4-1-1 dimensions motivates a model in which a complex scalar 
field in 4-1-1 dimensions, with the quantum numbers of the quark bilinears qu^Rj with fiavor 
labels i,j, fiuctuates about a background configuration related to the quark mass (the source 
of qq) and chiral condensate (the expectation value of qq). We will review this version of 
the hard-wall model in more detail in Sec. HIl As a demonstration of the pattern of chiral 
symmetry breaking in the model, the Gell-Mann-Oakes-Renner (GOR) relation for the pion 
mass was shown to be approximately satisfied jsl. Il3|. 
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where m^r and f^^ are the pion mass and decay constant calculated from the model, and rriq 
and a are parameters in the scalar field background playing the role of the quark mass and 
chiral condensate, respectively. The GOR relation is also satisfied in an SU(3) extension of 
the hard- wall model with an independent strange quark mass parameter 14 1. 

Classical calculations in the hard-wall model have reproduced a variety of QCD observ- 
ables with surprising accuracy, generally at the 10-15% level |5|, l6|. The hard- wall model fails 



at high energies 15|] where the Kaluza-Klein spectrum diverges from the Regge spectrum 



16|, a problem partially corrected in the soft- wall model 17|. More surprisingly, it was 
discovered that pion condensation in the hard-wall model has qualitatively different features 
from predictions of the chiral Lagrangian If the isospin chemical potential increases 

beyond a critical value, hadronic matter is expected to undergo a phase transition to a state 



in which a linear combination of the pion fields condenses |19l-l2l|. In chiral perturbation 
theory without unphysically large low-energy coefficients, the pion condensation transition 
is second order and approaches the Zel'dovich equation of state for stiff matter smoothly 
across the phase boundary [22]. The pion condensation transition has been studied using 



other approaches including lattice QCD 23|-|26|, with results in agreement with chiral per- 



turbation theory. However, in the hard-wall model the transition from the hadronic phase 
to the pion condensate phase was found to be first order, rapidly approaching the Zel'dovich 
equation of state across the phase boundary. The holographic model takes as input the 
pattern of chiral symmetry breaking, so disagreement with lowest-order chiral perturbation 
theory is surprising. 

Another more subtle puzzle related to chiral symmetry lies in the form of the GOR 
relation when the chiral condensate a is allowed to have a phase with respect to the quark 
mass. As deduced from a linear sigma model, the quark mass and condensate are in phase, a 
fact related to the existence of an anomalous chiral U(l) symmetry. However, beginning with 
the non-Abelian nonlinear sigma model we can ask what would happen if the parameter a 
corresponding to {qigR) were complex. In unpublished work, it was found that the resulting 
GOR relation as derived in the hard-wall model with complex a disagrees with the analogous 
prediction of the chiral Lagrangian. In particular, in the hard-wall model the GOR relation 



takes the form 
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mlfjog^^)=2m,{a-a*), (1.3) 

while the corresponding prediction based on the chiral Lagrangian and PCAC is the same 
as Eq. flL2p with a replaced by its real part. It is the goal of this paper to reconcile these 
discrepancies and restore consistency between holographic QCD and chiral perturbation 
theory. 

We will show that the incorrect structure of pion interactions in the hard-wall model is a 
result of the choice of boundary conditions imposed on the 5D fields. Care must be taken in 
order that the boundary conditions respect the symmetry breaking structure. The subtlety 



compared with other extra-dimensional models is that here the background of a bulk field 
not only spontaneously breaks the gauge invariance in the bulk, but a non-normalizable 
term in the background also explicitly breaks the gauge invariance. In the hard-wall model 
the bulk gauge invariance is responsible for the global chiral symmetry of the effective 4D 
theory, so it is important that the explicit and spontaneous breaking be correctly accounted 
for at the boundary. 

We will demonstrate that nonlinear boundary conditions on the bulk scalar field, or else 
nonlinear boundary terms in the action, may be consistently chosen so as to restore the 
proper pattern of chiral symmetry breaking in the hard-wall model. The nonlinear bound- 
ary dynamics we propose are an alternative to the description in Ref . j28| , which also leads to 
acceptable symmetry structure (and also accommodates a bulk Chern-Simons term absent 
in the present model). The nonlinear boundary conditions relate the 5D a-model scalars to 
products of pseudoscalars tt^'tt", (vr^vr")^, etc. In order to motivate these unusual boundary 
conditions and to demonstrate their relation to a Sturm-Liouville system, as required for 
consistency of the standard Kaluza-Klein decomposition of the fields and their interactions, 
we reparametrize the 5D fields by a nonlinear field redefinition. The reparametrization in- 
troduces a new surface term [i.e. a total derivative) involving the pions in the 5D action 
but replaces the nonlinear boundary conditions with ordinary linear boundary conditions 
consistent with the desired symmetry-breaking pattern. As opposed to the boundary con- 
ditions proposed to describe multitrace operators in the AdS/CFT correspondence 29|, l30| . 
the nonlinear boundary conditions in the hard-wall model arise away from the boundary of 
AdS, at the infrared boundary of the spacetime. 

The modifications of the hard-wall model as described in Ref. jsl required to restore the 
structure of chiral symmetry breaking have a number of phenomenological consequences. 
The GOR relation for the pion mass, Eq. (11. 2p . is correctly normalized only after the quark 
mass and chiral condensate are rescaled. This same rescaling is consistent with the AdS/CFT 
correspondence, and is a result of the modified boundary physics. The pion potential is 
modified so as to reconcile properties of the pion condensation transition with predictions of 
chiral perturbation theory. Most hadronic observables receive only small corrections which 
vanish in the chiral limit, so the hard- wall model remains surprisingly accurate in many of 
its predictions for low-energy QCD observables. 
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II. REVIEW OF THE HARD- WALL MODEL 

Following the conventions of Ref . jsl , the hard- wall model is defined by the 5D action 

S = l d'x^lTr [{Fi^lf + (Fiflr) + \DmX\' + , (2.1) 

where F^^) = F^^^''a^/2 and F^^) = F(^)V72 are field strengths for the 5D SU(2)lxSU(2)r 
gauge fields; a"" are the three Pauli matrices; M, G 0, . . . , 4 are Lorentz indices contracted 
with the AdSs metric from Eq. (11. ip : and X is a 2x2 matrix of complex scalar fields 
transforming in the bifundamental representation of SU(2)ixSU(2)/j. For the calculations 
in this paper we work with R= 1. 

The equations of motion have a solution with vanishing gauge fields and scalar field profile 

^o(^) = \{m,z + a/)l = ^1, (2.2) 

where U is the 2x2 identity matrix. The fields X have the quantum numbers of the scalar 
quark bilinears, which are the operators whose coefficients in the Lagrangian of the 3+1 
dimensional theory are quark masses. We approximate isospin as unbroken, so that the up 
and down quarks have equal mass. The term in the solution Eq. (12. 2p proportional to z 
is non-normalizable and is related by the AdS/CFT dictionary to the quark mass, which 
explicitly breaks the chiral symmetry; and the term in the solution proportional to z^ is 
normalizable and is related to the condensate {qiQii)-, which spontaneously breaks the chiral 

symmetry HQ. 

The non-normalizable mode in the scalar field background explicitly breaks a bulk gauge 
invariance, but the presence of this mode is equivalent to spontaneous breaking due to a 
heavy Higgs field localized on the ultraviolet boundary {z = e) in the decoupling limit. To 
see this, we write the Higgs doublet (0^, 0°) as a matrix. 




H=\ I , (2.3) 



which transforms in the bifundamental representation of the chiral symmetry. In this form, 
the up and down quark Yukawa couplings take the form £Yuk = — ATr {i/^gLg^ + /i.e.}. 
Replacing {—)qL(lR with the 5D field X, the localized Higgs field has boundary action 

5e = ^ (fxTi\^\d^H\^ -V{H) + ^{HX^ + XH^)^, (2.4) 



where V{H) is the Higgs potential exhibiting spontaneous symmetry breaking. A similar 



coupling appears in bosonic technicolor models |31|, |32|. The factor of in the last term 
ensures the proper scaling with the field X near the UV boundary. We impose Neumann 
boundary conditions on X in the ultraviolet, modified by the presence of the boundary term 
(12. 4p . Replacing the Higgs field H by its vacuum expectation value (H) = (0*^)1, chosen 
real, the equations of motion and boundary condition for X are given by: 

(\d,x] - \dX + ^X = 0, (2.5) 



yZ^ J Z^ 

d,X\^ = -\{H), (2.6) 

where □ = 'q'^^d^du- Greek indices will always refer to 3+1 dimensions, and capital Latin 
indices will refer to 4+1 dimensions. 

By identifying the diagonal quark mass Mg = —X{H), the boundary condition becomes 
S^Xl^ = Mg. Near the boundary as e — ?■ 0, the solution for X consistent with this boundary 
condition takes the form X ^ MgZ. Thus, the coupling of X to a Higgs field localized 
at the UV boundary gives rise to the appropriate non-normalizable background solution 
for X, which justifies the presence of the non-normalizable background and its AdS/CFT 
interpretation as the source for the operator qq. However, the overall normalization of Mg 
in terms of the physical quark mass depends on the normalization of the field X, and may 
also be modified by infrared dynamics as will be the case here. In Eq. (12.21) we have set 



Mg = mg/2, although other normal 
products of scalar quark bilinears 



izations better match QCD predictions for correlators of 



33|. 



The fluctuations of X, which contain scalars and pseudoscalars (pions), are typically 
decomposed as 

Xoid{x,z) = (^^{mgz + az^)l + S{x, z)^ U{x,z), (2.7) 

where S" is a Hermitian matrix of scalars and U = exp[i'n'"'{x, z)a°'] is unitary. Any matrix 
can be written as a product of a Hermitian and a unitary matrix, and any Hermitian matrix 
function of x and z can be written as the term in parentheses in Eq. (12.71) . so this ansatz is 
completely general up to a U(l) factor relevant for the chiral anomaly but which will not be 
discussed here. 

The scalars and pseudoscalars decouple at quadratic order in the action, so we temporarily 
limit our attention to fluctuations with S = 0. In order to simplify the discussion we also 



(2.8) 



temporarily decouple the vector fields by taking = 0. We will include the gauge couplings 
in Sec. |Vl but they are an added complication which is not necessary to understand the 
main conclusions. 

The lightest pion Kaluza-Klein mode, 7r(a;, z) = 7f'^(a;, z)a^ /2 = Tc"'{x)ilj{z)a°' /2, has action 

= J d^'x Tr I (df.Ud'^U^ - A^\zY nixf^ | . 

As explained in Ref. js] and will also be explained in Sec. IIIIl the pion wavefunction is flat 
with iIj{z) ^ 1 except near ^ = e, so integrating over z yields the effective 4D action for the 
pions, 

Seff = J d^x ^Tr ^^d^Ud'^U^ - Amlfc^^ , (2.9) 

where is determined by the equations of motion and boundary conditions, and from the 
kinetic term we identify the pion decay constant 



dzah' 



(J 



(2.10) 



as e — 7- 0. The expression fl2.10p for also follows from an AdS/CFT calculation of the 

[ I 

transverse part of the axial vector current-current correlator [18j] . 



III. CHIRAL SYMMETRY BREAKING IN THE HARD- WALL MODEL 

The structure of the pion effective action (12. 9p demonstrates the discrepancy between the 
hard-wall model as defined above and chiral perturbation theory. The pion mass term in 
Eq. (12. 9p does not include the higher-order pion interactions required for the chiral symmetry 
to be maintained while the quark mass matrix transforms in the bifundamental represen- 
tation under the chiral symmetry (like a Higgs spurion). In the chiral Lagrangian the pion 
mass term is proportional to Tr {Mgll^ + M^U) , which displays the proper pattern of explicit 
and spontaneous chiral symmetry breaking. Beyond quadratic order in the pion fields, the 
hard-wall model as described above disagrees with the chiral Lagrangian, leading to unusual 



pion phenomenology inconsistent with chiral perturbation theory. 



terms in the pion potential in this context was also noted in Ref. 



35|. 



'he absence of quartic 
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Restoration of the correct pattern of chiral symmetry breaking may be achieved by mod- 
ifying the boundary conditions in a nonhnear way which mixes the scalar modes S and 
products of pseudoscalars vf'^Tr", as we discuss below. It will be convenient to rescale the 
quark mass parameter rhg = —2mg, so that the background profile of the field X takes the 
form 

Xo{z) = -rugz + — ^ ^1. (3.1) 
We then consider a nonlinear redefinition of the 5D fields as follows: 

X{x, z) = -rrigz + {^z^ + S{x, z)^ U{x, z), (3.2) 

which is to be compared with Eq. (12.71) . We write U{x,z) = exp{iiT"'{x, z)a°'). Now the 
pseudoscalar fluctuations in U{x, z) multiply the term in the background responsible for the 
spontaneous breaking of the chiral symmetry, but not the term responsible for the explicit 
breaking. With boundary conditions S{x,e) = '7i°'{x,e) = S{x,Zm) = 0, and a Neumann 
condition on vr" at Zm, the scalar and pseudoscalar modes again decouple and the pion 
action takes the form 

S = J d'^xTif^^id^Ud'^U^ -d,Ud,U^)^ + y"ci^a;Tr|^(t/ + f/t)|^^}, (3.3) 
where the last term is an IR localized boundary term due to a total derivative in the action. 
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which also 



The field parametrization Eq. (13. 2p is an alternative to those of Refs. 
lead to an acceptable model, but with a nonlinear term at the UV boundary z = e rather 
than at the IR boundary z = Zm- The Kaluza-Klein modes are solutions to the linearized 
equations of motion, which follow from the quadratic part of the action: 

/2 3 /* 
d^x^id^n^d^n" - d.Tr^d.Tr") - / d'^x mga{nV)\^^, (3.4) 

where 7r'^(x, z) = tt°'{x)iI!{z). The linearized equation of motion for the pion wavefunction is 

d.iz^d,^) = -mlz^^. (3.5) 

The Neumann boundary condition in the IR is modified by the boundary term in the action, 
with the result 



9zHz)\.^,^ = —rT^izJ- (3-6) 



Note that the boundary conditions here are linear, and the nonlinear boundary terms in 
the action are treated as interactions. In this form, the pions are described by a standard 
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Sturm-Liouville system. The solutions are in terms of Bessel functions and the normahzable 
solution has the expansion 

1 ^ h higher order in rriT^z j . (3.7) 

If mj^Zm <^ 1 then ip^z) 1 in the entire interval e < z < Zm- Substituting the expansion of 
'ip{z) into the boundary condition Eq. (13.61) . we find to leading order in rriq, 

2 4 

ml^^ = 2mga. (3.8) 

Using Eq. (I2.10p . which is not affected by the field redefinition, Eq. (13. 8p is just the Gell- 
Mann-Oakes-Renner relation 

mlf^ = '^niga, (3.9) 
justifying the interpretation of rriq and a as the quark mass and chiral condensate, respec- 



33| . Note that the quark mass. 



tively, up to a simultaneous rescaling of and a as in Ref. 
and in particular the product m^cr, is rescaled from the old parameter rhq and even has a 
different sign. This rescaling is required in order to obtain the correct normalization in the 
GOR relation, but is also consistent with the AdS / CFT interpretation of rriq as the source 
for the operator {—)qLQR whose expectation value is the chiral condensate. The condensate 
is obtained by varying the action with respect to the source rriq. Because of the additional 
boundary term, which scales as rriq, the quark mass needs to be rescaled with respect to the 
chiral condensate as above. 

We now derive the 4D effective Lagrangian for the redefined pions. Let us first focus on 
the 2;-derivative piece. Ignoring the higher KK modes and writing U = exp{i'n'"'{x)ip{z) a""), 
we find 



a^z^ 



dz — i)'{zy 



tt''{x)ti\x), (3.10) 



as before. Integrating by parts and using the equations of motion and the boundary condition 
for %p we arrive at: 

Sd I d^'x - ^ 



dz{-mla'^z^)tlj{zy + 2mqa ^'K^{x)n^{x). (3.11) 

Using the flatness of the ip{z) profile, the expression for in Eq. f l2.10p . and the GOR 
relation, the above expression in brackets vanishes. Including the boundary term in Eq. (13. 3p . 
the approximate 4D effective Lagrangian is equivalent to the lowest-order chiral Lagrangian: 

S = y"rf^xTr|^9^f/9'^f/t + ^^([/ + f/t)|. (3.12) 



As a consequence, the properties of the pion condensate phase and other aspects of pion 
physics now agree with the predictions of chiral perturbation theory thanks to the modified 
boundary dynamics. 



IV. NONLINEAR BOUNDARY CONDITIONS 



The nonhnear reparametrization of the bulk field X in Eq. (13. 2p allows for independent 
linear boundary conditions on the scalar and pseudoscalar modes while maintaining the 
proper pattern of chiral symmetry breaking. In terms of the original decomposition of X as 
per Eq. (12. 7p . the boundary conditions required to maintain the pattern of chiral symmetry 
mix the scalar and pseudoscalar fields in a nonlinear way. To understand the structure of 
the nonlinear boundary conditions we can expand the two field decompositions, Eqs. (12. 7p 
and (1321): 



X 



-{niqZ + az ) + S 



-nigZ 



(l + 2i7r - 27?^ + . . . ) 



(4.1) 



where vr = 7r"cr'*/2, and similarly for vr. Equating the ant i- Hermit ian parts of the two 
descriptions gives, to quadratic order in the fields. 



TT 



7f (^1 + 



n 1 



nig ^ 




1 r 


oz^ ) 


— 


rhg ^ 


+ 


w 


oz^ ) 


az'^ y 


rhg \ 


- 


rhg 


oz^ I 




a'^z^ 



{Sn — Sn) + {nS — nS) 



Sn — Stt { 1 + 



Sn + TiS 



m„ 



oz^ 



+ . . 
- h.c. 



+ 



(4.2) 



Similarly, the Hermitian parts give. 



S = S — rUgZ TT^ 



az^ (n^ 



TX 



) + i S (tt — tt) — h.c. 



+ . 



S + nXaZ n 



z-2 



g^ ,. "2 



Sn — nS 



(4.3) 



These expressions have been left in terms of the old mass parameter rhg, which is equivalent 
to — 2mq, and terms higher order in rhg have been dropped. 

In the new decomposition of the field X, the boundary conditions consistent with the 
chiral symmetry breaking structure are: 

rh„ 



Six, 



S{x, z„ 



= 7r{x, e) 
10 



0, d^TT{x,z)\, 



az-i, 



-TT. 



(4.4) 



In terms of the original decomposition of the field X, the boundary conditions are, 

5 = ^ S = -mqz ^2 + . . . , (4.5) 



S^TT = \tX d^TX = ^71 + (4.6) 



rriq Snio 
— a^TX = — ^ 

where the ellipses include terms higher order in the fields and in rhg, and terms that vanish 
when traced over. 

For these boundary conditions to be physically acceptable, the boundary variation of the 
action must vanish. Expanding the action (12.11) with = about the background as in 
Eq. (12.71) . we obtain 

3 1 ^^'^^ 

H — g- ^^5* + S^^ > + constant. 

where v{z) = rfiqZ + az^ as in Eq. (12.21) . The terms in Eq. (14. 7p with z-derivatives lead to 
boundary terms in the variation of the action. Expanding to quadratic order in S and vr, 
we find for the boundary variation of 5, 



5S = / d^xTi 



2 ~ 1 ^ 
rdzS [nig + Saz^) 



6S - 



2(mgZ + az^Y a ~ 
— avr 



z^ 



5^ \ . (4.8) 



To leading order in ifiq and in the fields, using 

5S = —fhqZ (tT ^TT + (JtT TT ) + . . . (4.9) 

from Eq. (14. 5p . we find that the boundary variation 6S indeed vanishes. The cancellation 
of boundary variations in this nonlinear fashion is novel in the context of extra-dimensional 
models, though it is reminiscent of the mixed boundary conditions of certain Higgsless 
models 37| in which the contributions to the boundary variation of the action from different 
fields cancel one another. 

In Kaluza-Klein theories an effective description of the lightest modes is often derived by 
simply neglecting the heavier modes and integrating the action over the extra dimension. 
Indeed, that is how we derive the pion effective action in this paper. However, consis- 
tency of this approach relies on orthogonality and completeness relations dependent on the 
Sturm-Liouville structure of the equations of motion and boundary conditions. It is a math- 
ematical question which classes of systems of differential equations with nonlinear boundary 
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conditions satisfy the completeness and orthogonality theorems of Sturm-Liouville systems. 
In holographic QCD we have seen that there is a nonlinear field redefinition after which 
the boundary conditions are of the linear Sturm-Liouville form. This justifies the effective 
description obtained by keeping only the lightest modes. 

V. COUPLINGS TO VECTOR AND AXIAL- VECTOR FIELDS 

Having derived the chiral Lagrangian in the (75 — )■ limit, we turn to the case of nonzero 
5D gauge couplings with dynamical gauge bosons representing the vector and axial vector 
mesons. Including the gauge fields, the action takes the form of Eq. (12. ip . The field X 
transforms as a bifundamental under the gauge group SU{2)i^ x S'f/(2)R and we will use 
the gauge fixing condition = i?" = 0. We will also be working mainly with the linear 
combinations = (L^ — -RJi)/2, the axial vector field, and = (L^ + -R^)/2, the vector 
field. The normalization of these combinations by a factor of 2 (rather than ^/2) is so that 
their kinetic terms are canonically normalized given the unconventional normalization of the 
gauge fields in Eq. (12. ip . 

We parameterize the fluctuations of the field X as in Eq. (13. 2p . To leading order the 
scalars and pseudoscalars are decoupled, so we focus only on the pseudoscalars and set 
5 = for the present discussion. As in the previous section, the boundary condition on the 
pion is modified as in Eq. (13. 6p . 

We will determine the pion decay constant as in Refs. by the residue of the axial 

current two-point correlator at = 0. The AdS/CFT calculation of the correlator is 
performed by taking two functional derivatives, with respect to the source of the axial 
current operator, of the action evaluated on the classical solution to the linearized equation 
of motion for the transverse part of A'^. The resulting correlator is in terms of the bulk- 
to-boundary propagator for A^, which is a particular solution to the transverse-projected 
linearized equation of motion. This equation of motion for A^(q', z)±_ is 



If we have a solution to Eq. (15. ip of the form yl^(g, 2;) = A{q, z)AQ^{q), with boundary 
conditions d^A^q, z)\z^ = and A{q, e) = 1, then A{q, z) is identified as the bulk-to-boundary 
propagator and A^ Aq) is the source for the axial current. The AdS/CFT prediction for the 




(5.1) 
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pion decay constant is then, 

/I ^4^^ . (5.2) 

95 ^ z=e 

In order to study the pions we note that the pion fluctuations identified in the field X 
mix with the longitudinal part of the axial vector field — )■ S^^*^, which has the same 
quantum numbers. The pions will be identified as the lowest mode of this coupled system. 
Since we have in mind a Kaluza-Klein decomposition of the fields and since, for the purposes 
of deriving the low energy theory, we are only interested in the lowest mode, we will make 
the substitutions tt"-{x,z) — t- 7r"(x)^/'(z) and (f)"'{x,z) — )■ 7r"(a;)0(z). The linearized equations 
of motion for ip{z) and (f){z) are 

"^91 \z J (5.3) 
m^(v0 — az^^jj) = dz ((7z^dzi>) 



v(j) - az^ip = —dz ( -dz 



where the fields satisfy the boundary conditions dz4>{z)\z^ = 0(e) = ipi^) = 0, and Eq. fl3.6p . 



A. Approximate Analytic Results 

We can obtain an approximate solution to the equations of motion, Eqs. (15. 3p . in the 
chiral limit, in a similar fashion to Ref. j5|. We find that the approximate solutions near the 
e boundary are (j){z) = 1 — A{0, z) and ip{z) = 0, while those away from the e boundary are 
(f){z) = 1 — A{0, z) and ip{z) ^ 1. Plugging the first equation of Eq. (15.31) into the second, 
approximating v{z) ~ crz^, and integrating once we arrive at 

ml f 1 ^ , 1 



2 1 -dz(l) dz 

ag^ \z z 



az-'dz^p. (5.4) 



Now if we evaluate this expression on the IR boundary, using our approximate solution for 
(f){z) and recalling the boundary conditions, we find 

1 dzA{0,z) 
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2mqail){zm). (5.5) 



By utilizing Eq. (15. 2 p and the fact that ip^Zm) ~ 1, we have once again derived the GOR 
relation. 

The derivation of the chiral Lagrangian mass term is similar to that of previous sections. 
The only contributions come from the term proportional to dzUdzW^ and the boundary term 
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proportional to U + W . In particular we have, upon integration by parts, 



+ |rf^xTr{^([/ + t/t)|^ 



m J 



ln^ix)n'^ix) 

(5.6) 



If we make use of the equations of motion, Eq. (15. 3p . and substitute our approximate solu- 
tions for (f) and tp, we find that the two terms in the square brackets cancel, to first order in 
rrig. Thus, to this order in nig, we are left with the mass term of the chiral Lagrangian: 



54d3 y"rf4xTr|^^(f/ + f/t) 



(5.7) 



B. Numerical Results 



We will now present a numerical analysis of the equations of motion, Eqs. (15. 3p . We choose 
a value e = 10~^ MeV for the UV cutoff, and determine the location of the IR boundary 
to be Zm = 1/(323 MeV), by setting the rho mass to 776 MeV [5|. We take = 271 as in 
Ref. ^1, noting that the derivation of this assignment is unaffected by our new choice for 
the form of the field X{x,z), with different background. With the values rriq = 2.36 MeV 
for the quark mass and cr = (333 MeV)^ for the condensate, we have = 140 MeV for 
the pion mass and /^r = 92.0 MeV for the pion decay constant. The solutions for tp{z) and 
(j){z) are plotted in Fig. [H along with their approximate solutions, namely ip{z) = 1 and 
(j){z) = 1 — A{0,z). The functions ipi^) ^md (f){z) are normalized to obtain a canonically 
normalized kinetic term in the low energy theory. The plotted numerical solutions illustrate 
the extent to which the approximations of the previous section are valid. 

We would also like to understand numerically how robust the GOR relation is in this 
model, with respect to varying some of the parameters. In particular, fixing 171^^ = 140 MeV 
(by adjusting rrig for fixed values of cr) we would like to see what happens for different values 
of a. Varying the condensate by sampling a discrete number of points from (290 MeV)^ 
to (360 MeV)'^, the pion decay constant takes values ranging from 79 MeV to 102 MeV 
and the quark mass goes from 2.52 MeV down to 2.21 MeV. In Fig. [21 we plot the ratio 
(m^/^)/(2mgCr) for the above specified range of values for the condensate. If the GOR 
relation holds, the ratio should be approximately 1 and we can see that the plot shows good 
agreement over the entire region. 
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(a) (b) 

FIG. 1: Figures (a) and (b) display the numerical solutions to Eq. (j5.3p for tjj{z) and (j){z), respec- 
tively. In both figures there are two curves plotted: the blue curves are the numerical solutions and 
the red, dashed curves are the approximate solutions (^(z) = 1 in figure (a) and </> = [! — A{0,z)] 
in figure (b)). These plots were made with nig = 2.36 MeV and a = (333 MeV)^. 



290 300 310 320 330 340 350 360 



7^ 



FIG. 2: The red long-dashed line is a horizontal line at 1. The blue dots are the values of the 



function, {m^ f^) / {2mqa) . For this plot, the quark mass is adjusted to fix 
point. 



140 MeV at each 



VI. CONSEQUENCES FOR HOLOGRAPHIC QCD 



We check that the change of background does not adversely affect some standard predic- 
tions of the hard-wall model; that is, we will compare our predictions with those of "Model 
A" in Ref. js]. We might expect the modifications to be marginal because only the pion 
physics should be sensitive to changes in rrig, but we have also changed a sign. Aside from the 
substitution rhg — )■ —2mg, the derivations of the equations of motion for the gauge fields and 
of the vector current correlators are unaffected by the new form of the X-field. Since the ex- 
pressions for the related observables we compute are likewise unchanged, we refer the reader 
to Ref. [sl, for more detail. In the following, we summarize the methods for obtaining the 

15 



results. To calculate the mass of the ai, we solve Eq. (15. ip for the normalizable mode ipai{z) 
with boundary conditions ipa-^i^) = dzipai{zm) = 0. The wavefunction ipai is normahzed 
such that A^{x) has a canonical kinetic term, where A^{x,z) = ipai{z)A^{x) . As derivable 
from the two-point vector current correlator, we use the following expression to calculate 
the decay constants in terms of the profile in the extra dimension: F^/^ = \/ip"{0)/ g^. And 
finally, by looking at the terms cubic in the fields, coupling Vnir, we make a prediction for 
QpTrn- AH relevant cubic terms are the following: 



S D e 



abc 



(fx V,^ 



zgt, z-^ ^ ^ ^ ' 



(6.1) 



where we have used the equations of motion for the vector field to obtain this expression. 
In order to calculate the on-shell (7p,r7r from the effective 4D theory, we integrate out the 
extra dimension and identify gp^^.^ as the coefficient of e"^^9'^7r"(x)p|^(x)7r'^(x), where pj^ is the 
lowest mode in the vector field KK decomposition. Extracting this coefficient from Eq. ( 16. ip 
we find 



El 
P 

J TT 



dz ipp{z) 



L{d,<pf + -^{v<p-az'^y 



(6.2) 



_zg^ ^ 

where the 2;-integral of the expression in brackets is normalized to The results are 
presented in Table H] and we find that the new predictions have not changed significantly 
compared to "Model A" - they are still on the 10% level, with the exception of gpmr- 



Observable 


Measured (MeV) 


Model (MeV) 


Model A (MeV) [5] 




140 


140 


140 


nip 


776 


776 


776 




1230 


1370 


1360 




92.4 


92.0 


92.4 




345 


329 


329 




433 


493 


486 


QpTTTT 


6.03 


4.44 


4.48 



a. 



TABLE I: The predictions of the model as compared with "Model A" and experimental values 
The results are based on fitting to rn-Tr, /t^, and mp, leading to the parameter choice ruq = 2.36 MeV, 
a = (333 MeV)3, and Zm = 1/(323 MeV). 
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VII. CONCLUSIONS 



Holographic QCD models are surprisingly successful in their predictions of low-energy 
QCD observables. However, it was discovered in earlier work that the pion condensation 
transition in one version of the hard-wall model has qualitatively different features than the 
predictions of chiral perturbation theory and other approaches. We have shown that this 
disagreement is due to a boundary effect related to the explicit breaking of the gauged chiral 
symmetry by the non-normalizable background of a 5D scalar field. To restore agreement 
with the chiral Lagrangian we modified the boundary dynamics, either by introducing non- 
linear boundary conditions on the fields, or by performing a nonlinear field redefinition which 
induced an infrared boundary term in the action. The field redefinition allowed us to relate 
the system with nonlinear boundary conditions to a standard Sturm- Liouville system which 
manifestly maintains the proper symmetry structure, justifying the subsequent Kaluza-Klein 
decomposition of the fields. The chirally improved hard-wall model makes predictions for 
low-energy QCD observables that agree with the original model to within 1-2%. 

It would be useful to further explore the consequences of the modified boundary physics 
with regard to pion observables. It would also be interesting to find additional applications of 
nonlinear boundary conditions to extra-dimensional model building, for example to Higgsless 
models and holographic technicolor models. Finally, the necessity for nonlinear boundary 
dynamics in the hard-wall model provides motivation for further study of the mathematical 
problem of differential equations with nonlinear boundary conditions. In particular, it would 
be useful to classify those systems of equations and boundary conditions that can be related 
to a Sturm-Liouville system by a change of variables. 
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